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Abstract. We establish a representation of a class of solutions of 3d Navier-Stokes 
equations in M.^ using sums over rooted trees. We study the convergence properties of 
this series recovering in a simplified manner some results obtained recently by Sinai and 
' other known results for solutions in spaces of pseudo-measures introduced initially by 

^SJ , Le Jan and Sznitman. The series representation make sense also in the critical case 

^ ' where there exists global solutions for small initial data and it allows the study of their 

^ ' long-time or small-distance behavior. 

S: 

■ 1- Introduction 

. We consider the NS equation in the form 

> 



O 







vtik) = e-\''\\o{k) + 1 e-l'^l'^*-^) / dk'{k,Vsik-k'))PkVsik')ds (1) 



where Vt G C(]R'^;C^) is the Fourier transform of the velocity field, (■, ■) is the scalar 
■ product in and Pk : ^ is the projection on the directions orthogonal to the 

O ■ vector k G M^, i.e. Pka = a — {k,a)k\k\~^. Eq. ((H) will be studied in the spaces 

^ ! a G [2, 3) and where v G <I>(a, a) if i; G C{R^; C^) with k ■ v{k) = and 

I ■ \\v\\a = sup < OO. 



We will write a = 2 + e with e G [0, 1). 

The spaces a) are interesting because, in general, they contain solutions with infi- 
^ ^ ■ nite energy and enstrophy so classical results about existence and uniqueness of solutions 

^ . do not apply. 

^1 In a series of papers, Sinai [3131101, studies eq. (H)) in these spaces giving elementary 

proofs that in a) with a > 2, there is existence of unique local solutions and that 
these solutions survive for arbitrary large time if the initial condition is small enough. 

Moreover in the "critical" space $(2,2) there is existence and uniqueness of global 
solutions for small initial data. This latter global result was initially proven by Le Jan 
and Sznitman [II] using a probabilistic representation (under the name of stochastic 
cascades) and afterwards reproved by Cannone and Planchon jS] in a more standard 
functional- analytic fashion. 

The analysis of the equation in various function spaces similar to $(a,a) is sum- 
marized in the work of Bhattacharya et al. in p. A more recent review of the current 
status of the stochastic cascades approach has been recently given by Waymire in [T^. 

We are interested in explicit series representations for these (local or global) solutions. 
When a > 2 Sinai [9] proved that the local solution can be represented by a absolutely 
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convergent series and in ^Ul he analyzed this series with the aim of understanding better 
the growth of the various terms. A different series representation appears also in the book 
of Gallavotti [7j. 

Our contribution will be to prove yet another series representation for the solutions in 
a) including the critical case a = 2 which was left open by the analysis of Sinai. 
This series representation is indexed by rooted trees. 

Rooted trees appear naturally in the series expansion of solutions to ODEs. They pos- 
sess remarkable algebraic properties which were masterfully exploited by Butcher jl] to 
provide a general theory of Runge-Kutta (R-K) methods for numerical integration. After- 
wards rooted trees appeared also in the work of Connes and Kreimer [B' on the mathemat- 
ical structure of renormalization in quantum field theory. The work of Brouder [21 El gives 
a short overview of the algebraic properties of rooted trees and explore some connections 
between R-K methods and renormalization. 

These widespread applications of rooted trees were the initial motivations for this work. 
In the following we show that rooted trees are a natural language in which the known 
results (and some new ones) about the eq. (jT} in the spaces a) can be proven rather 
easily. Moreover the representation with a series indexed by rooted-trees can be controlled 
in a straightforward way and provide informations on the solutions themselves, like the 
behavior for large times or for large wave-vectors. 

The plan of the paper is as follows. In Sec. |21 we introduce rooted trees which will be 
used in Sec. E]to prove the series representation for the solutions of NS equation. Next 
in Sec. 13.11 we make some observations on the different nature of some classes of terms 
which contribute to the series and which appeared originally in the work of Sinai |lUj . 
In App. El we review briefly, for sake of completeness, the question of existence and 
uniqueness problem for the equation (jH) in the spaces $(a,a;). At the end, App. iBl 
collects some proofs. 

2. Trees 

A rooted tree is a graph with a special vertex called root such that there is a unique 
path from the root to any other vertex of the tree. Here some examples of rooted trees: 



We draw the root at the bottom with the tree growing upwards (as real trees). In a 
rooted tree the order of the branches at any vertex is ignored so the following two are 
representations of the same tree: 



Given k rooted trees ri, ■ • • , we define r = [ri, ■ ■ ■ , r^] as the tree obtained by attaching 
the k roots of ti, ■ ■ ■ , to a new vertex which will be the root of r. Any tree can be 
constructed using the simple tree • and the operation [■■■], e.g. 



: V V 



[•]=: 




etc. . . 
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On trees we will define various functions. Denote T the set of all rooted trees and let 
I ■ I : T — i> N the map which counts the number of vertices of the tree and which can be 
defined recursively as 

I • I = 1, • • • = 1 + kil H h |rfc| 

moreover we define the tree factorial 7 : T — N as 

7(») = 1, 7([n,---,^fc]) = |[n,---,i-fc]|7(n)---7(^fc)- 

Last, we define the symmetry factor cr : T — N: this is defined recursively as a{») = 
1 and cr([ri • • ■ r„]) = s(ri, . . . , r„)cr(ri) ■ ■ ■ cT(r„) where s(ri, . . . , r„) is the order of the 
permutation group of the (ordered) n-uple (ri, . . . , r„) G T". In the sequel we will only 
need to consider the subset BT C T which contains rooted trees with at most two sons for 
each vertex, this is due to the bilinear nature of the non-linear term in the NS equation. 



3. Series representation 
If we let Ct{k) = the eq. ((H) above takes the form 

Jo Jr3 |k — k |°|/c'|" 

for function c G C(M^; C'^) such that sup^ ||q||o < 00, {k, c{k)) = and c(0) = 0. 
For simplicity write the above equation in the abstract form 

Ct = StCo+ / St-sB{cs,Cs) ds. (3) 
Jo 

where 5* is a bounded semigroup and B is the symmetrized bilinear operator B{c, d) = 
(fii(c, d) + Bi{d,c))/2 with 

,{k,c{k-k'))Pkd{k') 



Bi{c,d) =i\k\'' / dk 



\k - k'Hk' 



Let V = {c G C(M^C3) : c(0) = 0, {k,c{k)) = and ||c||o < 00} and for any T > 
define the Banach space Wt = Cb{[0,T],V) endowed with the sup norm. Define the 
bilinear operator B : Wt ® Wt Wr as 



Bt{c,d)= / St^sB{cs,ds) ds. 
Jo 



Lemma 3.1. For any a > 2 and any T > 0, the operator B is well defined and there 
exists an increasing function N : IR+ IR+ such that 

\Bt{c,d)\<NMo\\d\\o (4) 

where Nt tends to zero as t for any a > 2. Moreover when a = 2 we have a uniform 
bound sup^yQ Nt < N^: < 00 independent of T . 

Proof. The proof can be found in the paper of Sinai (Hj and consists in a direct estimate 
of the integral. Some general considerations on the bilinear operator are summarized in 

App. m □ 
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Now define the operator : BT x V ^ Wt by recurrence as 

(P{»]h) = B{S.h,S.h) <i){\T]-h) = 2B{S.h,<i){T-h)) (5) 

and 

0([rir2];/i) = 2i3(0(ri;/i),0(r2;/i)) (6) 
for any h G V,r, Ti, r2 G i3T wliere we let {S.h)t = Sth. 

To find solutions of eq. in Wt witli initial condition /i G V we set up Picard iterations 

G Wt}„ as uf^ = and 4"^^^ = Sth + Bt{u^''\u^''^). 

Lemma 3.2. The functions m^"-* /iawe the representation 

where BTn C i3T zs i/ie sei of rooted trees for which the leaves are at distance at most n 
from the root and where we conventionally let BT = 0. 

Proof. It is clear that the formula holds for n = (the sum does not contain any terms). 
Assume it holds for any k < n and let us prove it for n + 1: 



S.h + B(S.h,S.h) +2 V -^B(S.h,(^(r-h)) 



r^BT.-. ^^^^ 



s.h + <j>{.;h)+ Yl ^<^(M;^) 



cr r 



S.h 



(T r 



since cr([rr]) = 2cr(r)^ and o"([r]) = cr(r). □ 
The norm convergence of the series 

u = S.h+J2 -^<^ir;h) (8) 

reBT ^^^^ 

in Wt implies convergence of the Picard iterates u^"^ to the element u G Wt which satisfy 
eq. (jHl) in [0, T] with initial condition uq = h. 

Define the following function 6 : BT M: 

d{.) = 2, d{[T]) = 1 + e(r), d{[T,r,]) = d{T') + 0(r2) 
and note that h v-^ (ptsi'^', h) is an homogeneous function of order 6{t). Always holds 

(|r| + l)/2 < ^(r) < |r| + 1 (9) 
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as easily proven by induction on |r|. 

The following control of the coefficients of the series (jH)) is the main result of this note. 
Theorem 3.3. For any e G [0, 1) the following estimate is true 

where the constants satisfy: 

A A 

for some constant A depending only on a. 

The proof of this theorem is reported in App. iBl 
Remark 3.4. The constants C-r can he chosen as follows: 

Cr = A\^\-i{T)-''^. 

Now we can prove the following result about existence and series representation of 
solutions of eq. ((H) in the spaces (a = 2 included). 

Corollary 3.5. The series (0) has the following properties: 

a) for e G (0,1) and fixed ||/i||o it converges in norm for small t^, and solve the 
problem |^ in Wt. ; 

b) for e G (0, 1) and fixed T it converges in norm in Wt for ||/i||o small enough; 

c) fore = (i.e. a = 2) and for \\h\\o small enough it converges in norm in Woo o-nd 
define a global solution of the problem ([^. 



Proof. Using Thm. 13.31 and Remark 13.41 we find that there exists a constant B such that 

(10) 



\ut{k) - [Sth]{k)\ < ^i?Ma(r)-S(r)-^/'e-l'=l'*/(l^l+^)tl^l^/2||;,l|e(r) 



n>l 



where Z„ is the number of rooted trees in BT with n vertices. The following recursive 
relations can be used to bound the Z„'s: 



Zi — 1, Zn+l < Zn + ^ ] Zn^Z.^ 



«2- 

n-i+n2=n 



From this relation it is not difficult to prove that Z„ grows at most exponentially, i.e. 
there exists a constant D such that 

Z„ < D"(n+ l)-3/2. (11) 

(see for example 0, Sec. 3). 

Next, by induction we can prove that 7(r) > 21"^!^^. This bound is optimal since it is 
saturated by the binary trees for which every path from the root to the leaves has the 



(12) 
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same length. Using this bound, eq. ^ and eq. (fTT|) in eq. (fTUj) we get 

\ut{k)\ < +^Z„5"e-l^l'*/("+^)r^/l/i||i"+')/' (l A 

n>l 

< \\h\\, + Y,{DBf/^r{n + ly^'^ht^^'^'^ (l A \\ht^^'^" 

n>l 

SO the series (jH)) converges in norm whenever the geometric series 

n>l 

converges. This gives directly a),&),c). Indeed note that for e = the dependence on t 
disappear in this last series. □ 

Now, lets come back to the original variables. It is clear that the function Vt{k) = 
\k\~'^Ut{k) satisfy eq. in [0,T] when the series defining u converges in Wt- Here we 
are interested in the behavior of the global solutions when a = 2: 

Corollary 3.6. In the case a = 2 and when \\h\\2 is sufficiently small the global solution 
V of eq. (OJj with initial condition h has the following two properties: 

a) for fixed k G ]R^\{0}, lim^^oo \vt{k) \ = 0; 

b) for fixed t > 0, there exists two constants C3, C4 such that \vt{k)\ < C'se"'^*'^'^ as 
\k\ 00. 

Proof. By the same bounds performed in Cor. 13.51 we see that the function Vt{k) = 
\k\~'^Ut{k) satisfy the inequality 

\vt{k)\ < e-l'=l'*|/i(A;)| + J2cnn + l)-=^/Vl^l'*/("+i)||/il|^"+')/', k ^ 

n>l 

for \\h\\2 small enough to guarantee the convergence of the series 

n>l 

Then fixed k G M'^\{0} we have limt^oo e"''^'^*''*-"''^"'^^ = for each n and we obtain that 
\vt{k) I — >• as t — (X). 

Next, we want to estimate the series at fixed t and for — >^ 00 by Laplace method. 
Write 

n>l n>l 



The exponent in the sum of the r.h.s has a maximum for n ~ \k\Vt/\/\ log(Ci|' ^"^''^^ 
and so, when — > cx) we have 



n>l 



where C4 = 2/ \/ \ log(Ci||/i||2^^)| and C3 is some finite constant. □ 
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3.1. Remarks on particular subsets of trees. In the series dH} different classes of trees 
give different contributions. We define simple trees the trees with at most one branch at 
each vertex, i.e. of the form [•■■[•]■■■]• Short trees are instead trees for which at each 
vertex we have two branches, each of which carries a fixed proportion (a or 1 — a) of 
the vertices. Of course this will not be possible in general, so we allow the proportion to 
oscillate around a in the interval [a — Aa, a + Aa], for some fixed < Aa < min(a, 1 — a). 
Since rooted trees does not distinguish between branches at a vertex, we take here the 
convention that the branches are ordered by the number of vertices in the corresponding 
subtree. With this convention we can consider, without loosing generality, values of 
a e (0,1/2). 

We will denote BTq the set of simple trees and BT^ the set of short trees corresponding 
to the proportion a. 

We have a first simple lemma: 

Lemma 3.7. For r G BTq we have 7(r) = |r|!. For any a G (0,1/2) there exists 
constants Di, D2, D^, D4 such that, for any t G BTa we have 

^skr^M"' < 7(r) < Di|r|-iDf I. 

Proof. The proof of the first claim is trivial. For the second, note that we can choose the 
constants Di, D2, D3, D4 such that the inequalities are true for all the trees r G BTa with 
|r| < n for some fixed n and moreover they satisfy 



<1 and . . > 1. 



(a - Aa)(l - a- Aa- n'^) ~ (I - a + Aa){a + Aa) 

Then we proceed by induction on n>n. Assume the inequality is true for trees with 
|r| < and observe that, for r G BTa-, \t\ = n we have r = [T1T2] with 



and 



Then 



(a - Aa) < < (a + Aa) 
\t\ 



[1-a-Aa- \t\-^) < ^ < (1 - a + Aa 



|ri||r2| (a — Aa)(l — a — Aa — |r| ^) |r| 

~ (a — Aa)(l — a — Aa — n^-'-) |r| ~ |r| 

and similarly we obtain 7(r) > D^D^J^ItI^-^ , proving the claim. □ 

This different behavior of the two classes of trees is responsible for different convergence 
properties of the sum (jH} when restricted to simple or short trees. 
Define 

^-^ cr r 

then as consequence of Lemma f3. 71 and Theorem 13.31 we have the following result: 
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Corollary 3.8. For e > 0, the series Wt converges in V for every t and every initial 
condition /i G V and 

n=l ^ 

Proof. The estimates on the series are similar to those in Cor. but now the coefficient 
7(r) = |r|! goes to infinity fast enough to guarantee the convergence of the series for any 
time. □ 

In ^j, Sinai studied different classes of contributions to his series representation of 
NS. He calls the various contributions diagrams and then introduces short and simple 
diagrams which are analogous to short and simple trees (even if diagrams does not exactly 
corresponds to our trees). He then shows that the contribution of the simple diagrams 
cannot cause the divergence of the overall series. Corollary IH.8I is the analogous of this 
result in our setting. 

For short trees the function 7 behaves exponentially with the size of the tree and this is 
not enough to make the series restricted to short trees converge for arbitrary time (when 
e > 0). A similar phenomenon is observed in JU] for short diagrams. 
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Appendix A. Remarks on the spaces $(a,cj) 

Without going in detailed proofs we would like to note some remarks about the natural 
functional spaces in which solutions of eq. (P) live. Following Sinai |S| define the space 
^{a,u!) {a,uj > 0) as the space of continuous functions f : — > such that 

ll'^IU.a; = sup 1p{k)~'^\v{k)\ < OO 

where ip{k) = \k\~°' for \k\ < 1, ip{k) = \k\~'^ for \k\ > 1. Functions in this space can be 
bounded above by for small k and by \k\^'^ for large k. Consider the the bilinear 

integral operator 

8{v,v){t,k) = t [ e-l^l'(*-') /" dk'{k,Vs{k-k'))PkVs{k')ds (13) 



appearing in the r.h.s of eq. (^. For this operator we have the bound 

mv,v)it,k)\ < sup \\vs\\ljk\-\l-e-\'\'')lik) (14) 



0<s<t 



where I{k) = Jj^g dk'ipik — k')ip{k). For k the integral I{k) converges when u > 3/2 
and a < 3 and we prove easily that, when < 1 

for a < 3/2 
I{k) <Ca(\\og\k\\ hi a = 3/2 




while when > 1 

Hk) < a 



13-20.1 



for 3/2 < a < 3 



for 3/2 < < 3 



\k\-^ for > 3 



This behavior translates in the following estimates for the r.h.s. of eq. ()14j) . So when 
\k\ < 1: 

|^|-a^(i-a)/2 for a < 1 
|A;|-^(1 - e-l'^l'*)/(A;) < Ca { |A;|-"(1 A for 1 < a < 2 



|^|-aU|(a-2)/2 for 2 < « < 3 



and when Ifcl > 1 



^|2-2^(1 -e-l'^l'*) for3/2<cj<2 
\k\~\l - e~^''^"')I{k) < { \k\-'^{l A for 2 < < 3 

|A;|-^(1 A |t|)i/2 for3<cj 
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These bounds imply that B maps C([0, T], u;)) in itself whenever > 2 for any 
a G [0,3) and in this case the norm of i3 is given by 

Nt = snpsnp[tP{k)-^\k\-\l - e-l^l'*)/(A;)] 

t<T k 



and become small with T allowing a direct proof of existence and uniqueness of solutions 
to eq. (0) for small time. Moreover when a G [1, 2] the norm A^^- is uniformly bounded in 
T and this implies existence and uniqueness of global solutions with small enough initial 
condition. 

Note moreover that the same bounds are true on the torus (only wave- vectors |/c| > 1 
are important in this case) and that they always imply uniform control in time of the 
norm A^^ for any u > 2. In this case we have the existence and uniqueness of global 
solutions with small initial conditions whose decay at infinity is not worse than 
Details can be found in |8j. 



Appendix B. Proofs 

B.l. Theorem 



Proof. We will prove the statement by induction on |r| = n. Let us assume that the 
estimate is true for any tree r' with |r'| < n and let us prove it for trees r with |r| = n. 
Consider the case r = [t^t'^] with \t^\ = p, |r^| = q: 

t 

2 



i[Tir2];hm\<2 / e-l'^l (*-«)|5(0„(ri;/i),0„(r2;/i))|rfM 







<2 / rf^ze-l'^l (*-")|A;r+i / dk 







J \ Mri;h){k-k')\\Ur2;h){k')\ 
\k - k'Hk'l"' 



< 2C,^CAhC^ f due-\'\'^'--^u^'^^\^'\^\^'\^\k\^+' 

J s 



dk' 



^-\k-k'\-^u/{p+l)-\k'\^u/{q+l) 
„l\a 



\k - k'\°'\k\ 

The exponent in the integral has a maximum as a function of k' and 



\k — k'\'^u \k'\'^u Ifcpu 
> 



p+1 g + 1 p + g + 2 
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for any fc' G M^. So 

\M[rM,h)ik)\<2C,^C,2\\hC'> 

. g-|fc|2t/{P+q+2) r^yg-|fcP{t-«)(p+g+l)/(p+g+2)(^^^e/2(|rl| + |r2|)|^|a+l 

Jo 

dk'-. 



\k - A;'|°|A;'|" 

^yg-|fcP(i-«){P+'?+l)/{P+'?+2)(^^)e/2{|ri| + |r2|)|^|4-a 

where, if e is a unit vector in M.^ we let 

A' = 2 [ dk'- ,, , < cx). 

Consider the term on this last line of eq. (fT^ : 

^yg-|fcP(<-«)(P+g+l)/(P+'?+2)(^)e/2(|Ti| + |r2|)|^|4-a 
= ^(|rl| + |r2| + l)e/2 /"^ | ^ 1 4-ag-|fc|2 (l-«)(p+g+l)/{p+g+2) ^e/2{|Tl | + |r2|) 



(15) 



with A; = {t — sY^'^k. Let a = (p + g + l)/(p + g + 2). We have the following bound, proved 
in lemma IRTI below: 

1 - / ^ \ -1 

g-|fcp(i-.)a^b+g)e/2^^ < + {p + q)e/2 

Gathering all together we get 

mnr,];hm\ < Au.aHi/^iir^e-i^iv(ki+i),ki^/2 i^i'^' 

a\k\^ + {p + q)e/2 

When e G (0, 1) we have 

fc6K3 a|fcr + (p + g)5/2 
where iiT is a constant not depending on e, so in this case 

When £ = we have instead the bound 



r 



e/2 



Proving the claim in this case. The other cases can be treated similarly. □ 
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Lemma B.l. 



Proof. Easy: 



Jo 

Jo Jo Jo 

Jo 
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